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摘 要 : 为 了 揭示 河西 内 陆 河 浅 山 区 小 流域 的 蒸 散 发 及 干旱 特征 ,以 敦煌 西 土 沟 流域 为 例 ,利用 常规 气象 数据 , 通 
过 Penman-Monteith 公 式 计算 西 土 沟 流 域 8 个 观测 站 点 的 洪 在 蒸 散发 ,再 此 基础 上 利用 双 作物 系数 法 计算 出 各 个 观 


测 站 点 的 作物 系数 ,由 此 计算 出 实际 燕 散 发 ,并 计算 了 该 流域 的 SPEI( 标 准 化 降水 莹 散发 指数 )。 结 果 表 明 :(1) 西 
土 沟 流 域 年 潜在 蒸 散发 量 和 实际 蒸 散发 量 分 别 为 978 mm .258 mm; 上 游 到 下 游 随 海拔 下 降 ,潜在 蒸 散发 和 实际 蒸 


散发 随 之 呈现 不 同 程度 的 上 升 趋势 ;(2) 该 流域 下 游 干 旱 持 续 时 间 较 长 ,整个 生长 季 都 会 出 现 不 同 程度 的 干旱 , 秋 
季 是 干旱 发 生 的 主要 季节 ;(3) 流域 上 游 的 降水 量 大 于 实际 蒸 散发 量 ,中 下 游 实 际 蒸 散发 量 大 于 降水 量 ;(4) 利用 
NDVI 值 计算 的 作物 系数 玉 值 在 实际 蒸 散发 的 估算 中 有 较 好 的 适用 性 ,在 内 陆 河 浅 山区 小 流域 的 干旱 评估 中 ,干旱 
指数 SPEI FE SPI( 标 准 化 降水 指数 ) 和 PDSI( 帕 尔 默 干旱 指数 ) 指 数 更 有 优势 。 本 研究 计算 出 了 西 土 沟 流域 的 潜在 


和 实际 蒸 散发 ,并 评估 了 该 流域 的 月 尺度 及 季节 尺度 的 干旱 特征 ,为 该 流域 特别 是 中 下 游 的 生产 生活 用 水 提供 指 


导 , 并 得 到 适用 于 整个 河西 内 陆 河 浅 山区 无 资料 河流 的 普 适 性 规律 。 
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河西 内 陆 河 浅 山区 面积 大 、 范 围 广 ,其 浅 山区 
面积 约 占 总 山区 面积 的 11.2%。 河 西 内 陆 河 浅 山区 


究 , 对 农业 用 水 管理 和 区 域 水 资源 优化 配置 具有 重 
要 意义 ”。 实 际 蒸 散发 是 指 在 一 定 的 气象 条 件 下 ， 


出 山地 表 水 资源 量 为 0.2197x10: mi, 对 整个 河西 地 
区 有 着 重要 的 意义 。 该 地 区 除了 有 实测 气象 水文 
数据 的 较 大 型 的 河流 外 ,其 他 河流 无 实测 的 气象 、 
水 文 等 资料 ,这 些 河流 虽 水 量 不 大 ,但 却 有 独立 的 
供水 意义 ,下 游 孕 育 着 不 同 的 绿洲 ,对 当地 的 自然 
环境 及 人 类 生产 生活 有 重要 影响 。 因 此 ,本 研究 选 
取 典 型 的 浅 山 区 无 资料 流域 一 一 敦煌 西 土 沟 流域 ， 
对 该 流域 的 降水 、 蔡 散发 和 干旱 特性 进行 研究 , 揭 
示 该 流域 的 水 量 平衡 ,以 期 得 到 适用 于 整个 河西 内 
陆 河 浅 山区 无 资料 河流 的 普遍 性 规律 。 

蔡 散 发 是 地 表 水 量 和 热量 平衡 的 重要 参量 ,也 
是 水 循环 和 能 量 循环 的 重要 组 成 部 分 ,影响 着 区 域 
75 BER AE AS .气候 和 农业 生产 六 。 洪 在 蒸 散发 
是 指 在 一 定 气 象 条 件 下 水 分 供应 不 受 限 制 时 , 陆 面 
可 能 达到 的 最 大 蒸 散发 量 。 潜 在 蒸 散发 作为 农田 
灌溉 管理 .作物 需 水 量 估算、 稀缺 资料 地 区 水 量 平 
衔 等 研究 中 的 重要 参量 ,其 时 空 分布 及 变化 规律 研 


收 稿 日 期 : 2019- 12 -23; 修订 日 期 : 2020-05-14 


基金 项 目 :国家 重点 研发 计划 项 目 (2016YFC0400908); 内 蒙古 自治 


下 垫 面 受 水 分 供应 的 影响 而 实际 所 能 达到 的 蒸 散 
发 量 ”。 主 要 包括 两 个 方面 ,一 是 发 生 在 土壤 表面 
或 水 面 的 蒸发 ,二 是 发 生 在 植被 表面 的 蒸腾 ,实际 
蒸 散发 的 计算 通常 要 包括 这 两 个 方面 以 反映 当地 
的 实际 蒸 散 发 情况 。 

关于 薰 散发 的 研究 较 多 , 主要 集中 在 以 下 儿 个 
方面 :Q@ 潜在 蒸 散 发 (ETu) 时 空 演变 研究 ,比如 利用 
Penman-Monteith 公式 2 | Hargreaves 模型 "9 priestley- 
taylor #7” Blaney-Criddle 模型 5 、Rohwer 模 型 ”等 
方法 计算 潜在 蒸 散发 ;@) 实际 燕 散 发 的 估算 , 比如 
利用 FAO 56 推 荐 的 作物 系数 法 (K.) ,实际 车 散发 
ET=EToxK.*"”" 利用 涡 度 相关 法 “或 蒸 渗 仪 法 估 
算 ” 基 于 经 验 统计 或 水 量 平衡 方法 的 遥感 模型 "; 
© 作物 蒸 散 发 的 驱动 机 制 ,气候 变化 和 人 类 活动 对 
莹 散发 的 影响 “” ;中 莹 散发 计算 方法 的 区 域 适 用 
性 评价 2 。 在 我 国 干旱 地 区 ,由 于 土壤 水 分 的 限 
制 和 植被 生长 状况 的 差别 ,实际 蒸 散发 显著 小 于 潜 
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在 蒸 散发 ,所 以 在 河西 内 陆 河 流域 ,如 果 以 潜在 蒸 
散发 量 作 为 流域 水 量 的 支出 要 素 ,会 产生 较 大 的 误 
差 ,无 法 指导 人 类 的 生产 生活 取 用 水 ,因此 ,对 当地 
实际 蒸 散 发 的 估算 就 显 的 尤为 重要 。 但 是 由 于 实 
际 蒸 散发 的 估算 难度 较 高 ,所 以 ,以 往 的 研究 中 关 
于 洪 在 蒸 散发 的 研究 居多 ,关于 实际 蒸 散发 的 研究 
相对 较 少 ,在 流域 尺度 的 实际 蒸 散 发 估算 研究 更 
少 。 因 此 ,本 研究 希望 可 以 以 河西 内 陆 河 浅 山区 小 
流域 为 研究 对 象 ,通过 对 该 流域 上 游 到 下 游 不 同 海 
IE ,不同 流 域 位 置 的 潜在 和 实际 蒸 散发 进行 估算 ， 
得 到 适用 于 整个 河西 内 陆 河 浅 山区 无 资料 河流 的 
普 适 性 蒸 散 发 规律 。 

蒸 散 发 是 水 文 气象 研究 中 的 关键 变量 ,其 在 不 
同 海拔 和 不 同 植被 类 型 区 域 对 气象 要 素 的 响应 特 
征 不 尽 相 同 。Penman-Monteith 模型 是 世界 粮农 组 
织 (FAO)1998 年 推荐 使 用 的 潜在 蒸 散发 估算 方法 ， 
该 模型 是 目前 公认 的 无 论 在 湿润 还 是 干旱 地 区 计 
算 潜 在 蒸 散 量 精 度 较 高 的 方法 之 一 。 所 以 ,本文 在 
总 结 前 人 研究 结果 的 基础 上 ,首先 使 用 Penman- 
Monteith 公 式 计 算 潜 在 蒸 散发 ETv, 后 采用 文献 [24 ] 
和 文献 [25] 的 研究 中 得 到 的 作物 系数 到 与 NDVI 的 
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图 例 
(a) 我 国 西北 部 分 地 区 
(b) 西 土 沟 流域 

(c) 下 游标 柳 林 气象 站 
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关系 式 来 计算 不 同 作物 类 型 ,不 同 生长 期 和 非 生长 
期 的 作物 系数 K., 并 将 土壤 水 分 胁迫 因子 K. 考 虑 其 
中 ,最 终 得 到 不 同 海拔 区 域 流域 上 中、 下 游 实际 燕 
散发 。 此 外 ,将 基于 降水 和 潜在 蒸 散发 的 干旱 特征 
指数 考虑 其 中 ,以 期 对 整个 无 资料 的 河西 内 陆 河 浅 
山区 小 流域 的 蒸 散 发 和 干旱 特征 研究 提供 参考 , 进 
一 步 为 此 类 型 的 小 流域 的 水 文 气象 研究 提供 基础 。 


1 研究 区 概况 


西 土 沟 流域 发 源 于 阿尔 金山 北 项 ,最 高 海拔 
5536 m, 由 红柳 沟 、 阿 克 塞 沟 、 黑 沟 等 数 条 常年 性 流 
水 沟谷 组 成 (图 1) ,上游 以 中 高 山 山地 为 主要 地 貌 
单元 ,海拔 4100 m 以 上 发 育 现代 冰川 。 中 游 以 洪 积 
扇 为 主要 的 地 貌 类 型 ,下游 地 貌 类 型 相对 复杂 , 包 
括 沙漠 KEE .盐碱地 rc A RT 
域 面积 约 1125 km? ,年 均 地 表 水 径流 量 0.1104 10° 
m。 植 被 呈 地 带 性 分 布 ,在 流域 上 游 ,从 高 海拔 到 
低 海拔 分 别 为 亚 冰雪 带 ,高山 灌 丛 .草原 荒漠 带 、 山 
地 草原 带 。 流 域 中 游 主要 流 经 阿 克 塞 盆地 ,植被 以 
小 灌木 为 主 ,自然 景观 以 荒漠 戈壁 为 主 , 植 被 覆盖 


注 :(a) 我 国 西北 部 分 地 区 及 西 土 沟 流域 的 相对 位 置 ;(b) 西 土 沟 流域 及 流域 内 气象 站 的 分 布 位 置 ; 
(@) 流 域 下 游 架 设 在 树 柳 林 中 的 气象 站 ;(d) 流 域 上 游 架 设 在 青 崖 子 山 上 的 气象 站 ,海拔 3370 m, 
图 1 西 土 沟 流 域 地 理 位 置 和 流域 内 气象 站 点 分 布 示 意图 


Fig. 1 Geographical position of the Xitugou watershed and the distribution of meteorological stations in the watershed 
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度 在 20%~40%。 流 域 下 游 植被 以 荒漠 灌 丛 荒漠 草 
地 农田 为 主 。 本 区 地 处 戈壁 沙漠 区 ,地 势 比较 平 
坦 ,干旱 少雨 ,植被 覆盖 度 低 , 属 典型 的 大 陆 性 气 
候 。 其 特点 是 :夏季 炎热 .冬季 寒冷 .昼夜 温差 大 、 
降雨 少 .蒸发 量 大 空气 干 燥 、 积 温 高 .日照 长 SE 
低 。 区 内 年 平均 日 照 时 数 为 3240 h, 年 平均 气温 
9.4% ,最 冷 (1 月 ) 为 -9.2 % ,年 气温 差 34.1 % ,年 平 
均 无 霜 期 145 d ,绝对 无 霜 期 100 d。 多 年 平均 降水 
量 40.1 mm, 降 水 量 年 内 变化 较 大 , 主要 集中 在 6 一 8 
月 ,其 降水 量 约 占 全 年 降水 量 的 75% ,多 以 暴雨 出 
现 。 蒜 发 严 测 量 的 多 年 平均 蒸发 量 高 达 2417.19 
mm, 约 为 降水 量 的 60 倍 2 。 


2 数据 与 方法 


2.1 数据 来 源 

在 西 土 沟 流域 上 、 中 下游 共 架 设 8 套 小 型 自动 
气象 站 (型 号 :WPH1/HOBO-U30-CSM) ,在 每 个 气 
象 站 下 额外 架设 了 8 套 土 壤 水 分 观测 探头 ,水 分 探 
头 埋 深 分 别 是 10 cm、20 cm、40 cm、60 cm, 80 cm. 
站 点 名 称 分 别 为 敦煌 宫 (DHG) FEK 6 (BST) 青岩 
子 沟 (QYZG) 青 崖 子 山 (QYZ3) PEI AK (CLL) 、 阿 
TEIE (AKS) , ILRI) OLH) 、 青 岸 子 项 (QYZ4)8 个 观 
测 站 点 (图 1), 以 下 使 用 站 点 缩写 ,基本 信息 如 表 1 
所 示 。 气 象 数 据 来 自 于 在 上 述 观 测 点 架设 的 小 型 
气象 站 的 日 最 高 气温 .日 最 低 气温 .日 均 温 、 相 对 湿 
度 . 2 处 风速 日照 时 数 .气压 ;土壤 水 分 数据 是 由 
每 个 气象 站 额外 铺设 的 5 层 水 分 探头 测量 得 出 。 气 
象 和 土壤 水 分 数据 均 以 30 min 为 观测 间隔 ,将 其 经 
过 异常 值 吻 除 后 处 理 为 日 值 数据 ,数据 时 长 为 2018 
年 5 月 1 日 至 2019 年 4 月 30 日 的 数据 。NDVI 数 据 
(MOD13A1 数据 集 ) Æ A NASA (https://ladsweb. 
modaps.eosdis.nasa.gov/) ,数据 空间 分 辨 率 为 250 m, 
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时 间 分 辨 率 为 8 do 

2.2 研究 方法 

2.2.1 Penman-Monteith 公 式 ” 世 界 粮农 组 织 (FAO) 
在 1998 年 给 出 的 修正 Penman-monteith 方 程 是 目前 
被 广泛 使 用 的 车 散发 估算 方法 之 一 ,其 参考 下 热 面 
为 高 度 0.12 m 的 供水 条 件 不 受 限 制 的 草地 ,计算 
公式 为 : 


0.408A (R, - G) Ey U 
a A+y(l +0.34u,) 

SU LET) AR EZR AL (mm: dd!) ;6 为 土壤 热 通 量 
(MJ-m?+d"); 7242 m fea BEAD AY H IJAH CC); UA 
2 m 高 度 处 风速 (ms ');e. 为 饱和 水 汽 压 (kPa);e, 为 
实际 水 汽 压 (kPa) ;ee 为 饱和 水 汽 压 差 (kPa) ;A 为 
水 汽 压 曲线 斜率 (kPa°C7);y 为 干 湿 表 常数 (kPa*°C7); 
RASS (MJ +m? +d") ,是 净 短 波 辐 射 (R,) 和 净 长 
波 辐 射 (R,) 的 差 值 ( 式 2) ,由 于 架设 的 小 型 气象 站 
没有 观测 日 照 时 数 , 所 以 本 研究 中 采用 Hargreaves 
的 辐射 公式 进行 太阳 短波 辐射 RR 的 估算 ( 式 4): 
R =R -R (2) 

R,=(1-a)R, (3) 


R, = Kgs T max 7 Tmin) Ra (4) 
T ax. s F T iia j 
R= on -0.14 Je,) 


R, 
1.38 R 0.35 


so 


(5) 


R =(a,+b)R, (6) 
SUH: R WPK (MJ em? d); Ra D 
I (MJ-m?+d") sR KAUE TK BASE M-m’. d); 
7 和 Ts 分别 是 日 最 高 气温 和 日 最 低 气 温 。 扩 ,为 
经 验 调整 系数 , 取 值 0.16-~0.19 ,单位 是 "C" ,在 内 陆 
地 区 ,Rs=0.16 ,沿海 地 区 Kn.~0.19。R. 表 示 太 阳 短 
波 辐射 (MJ'm”.d); 及 表示 晴空 太阳 辐射 (MJ'm”， 


表 1 观测 站 点 信息 


Tab.1 Basic information of observation stations 


名 称 纬度 (N) 经 度 (E) 海拔 /m 下 垫 面 
JLŽEWI (JLH) 39°55'37.30" 94°0'20.98" 1236 沙漠 
敦煌 宫 (DHG ) 39°55'45.94" 94°1'29.62”" 1249 葡萄 地 
拜 水 台 (BST) 39°54'58.16" 94°1'55.53" 1273 荒漠 草地 
标 柳 林 (CLL ) 39°54'33.40" 94°2'45.81" 1274 人 工 树 柳 林 
阿 克 塞 (AKS) 39°38'5.75" 94°1'13.69" 1686 HEAR A 
FE FY (QYZG) 39°21'47" 94°1'11.04" 3100 浅 山 区 草地 
青 崖 子 3(QYZ3 ) 39°2116.13" 94°0'54.18" 3370 高 山 灌 从 草地 
青 崖 子 4(QYZ4) 39°18'15.76" 93°57'37.92" 4100 寒 漠 带 
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dd ) ;RJ/R,, 为 相对 短波 辐射 (其 值 <1.0); 由 于 没有 实 
测 的 太阳 辐射 资料 ,a, 和 5b. 值 采用 FAO 56 中 建议 的 
取 值 ,a.=0.25 ,45.=0.50。R, 的 计算 方法 如 式 7 所 示 : 


aa 


R,=——G,,d,[w, sin(g)sin(6) + 


(7) 
cos on (6)sin (w,)| 

式 中 :6 为 太阳 磁 偏 角 ;q 为 地 理 纬度 ;G. 为 太阳 常数 = 
0.0820 MJ.m ”min ;id 为 日 地 相对 距离 的 倒数 ;mw、 
为 太阳 日 落 时 角 。 详 细 的 计算 方法 可 查阅 FAO 56%, 
2.2.2 作物 系数 K。 双 作 物 系数 法 ”: 关 于 作物 系 
数 玉 的 双 作 物 系数 法 ,是 由 Wright 在 1982 年 提 
出 ,这 种 方法 将 作物 系数 分 为 植物 蒸腾 系数 (Ku) 和 
THR RAR BUK.) o 


K =KK, +K, (8) 
(NDVJ),,,. -NDVD vst 
= 0.54 
SORUN ， - oe T | 
K=BA=j,) (10) 
f.=1.3514x NDVI - 0.2811 (11) 


式 中 : 开 为 作物 系数 ;KK 为 基础 作物 系数 ,表征 作物 
蒸腾 的 影响 程度 ;及 为 土壤 莱 发 系数 ,表征 土壤 蒸发 

的 影响 程度 ;及 为 水 分 胁迫 系数 ,表征 根 区 土壤 含水 
率 不 足 时 对 作物 蔡 腾 的 影响 ,K. 的 计算 是 根据 土壤 
含水 量 和 FAO 56 规 定 的 各 类 土壤 的 田间 持 水 量 和 
凋 雁 系数, 通过 分 段 函数 计算 得 出 。NDVI. 和 ND- 
VIww 分 别 代表 裸 土 和 植被 密集 处 的 植被 指数 ;NDVI 


ChinaXiv 合 作 期 刊 


Ww F 


为 植被 指数 ;f 为 植被 覆盖 土壤 表面 的 有 效 面积 比 ;B 
为 经 验 系 数 ,取决 于 作物 生长 初期 和 中 期 的 K. 值 3。 
2.2.3 标准 化 降水 蒸 散 指数 SPEI ”SPEI 是 在 SPI 的 
基础 上 ,以 水 循环 为 理论 基础 ,考虑 变化 环境 下 降 
水 和 潜在 茸 散 发 的 干旱 指数 。SPEI 计 算 过 程 是 先 
利用 改进 的 Penman-Monteith 方法 计算 潜在 蒸 散 发 
量 , 之 后 用 月 降水 量 与 月 潜在 蒸 散发 量 的 差 值 作为 
SPEI 的 输入 数据 ,最 后 以 3 参数 (尺度 参数 .形状 参 
数 和 Origin 参数 ) 的 log-logistic 分 布 的 频率 值 计 算得 
出 SPEI 值 所 ,SPEI 指 数 的 干 湿 等 级 划分 所 如 表 2 
所 示 。 


3 cote Wc, W? 

1+dW+d, W +d W? 
W=J-2InP ,P<0.5 A 
Tepe. cyte W +e, W? 

1+dW+d, W’ +d,W ) 

W= |-2In(1-P), P>0.5 
SUF : co, c1, c2, di, da, ds YN FF BL, co=2.515517, ci= 
0.802853, c:=0.010328, di=1.432788, d=0.189269, 
ds=0.001308 ; Py BRIER ; WEE F P WI PR 


3 ARAM 


3.1 流域 潜在 蒸 散发 和 实际 蒸 散 发 的 日 变化 
由 图 2 可 知 ,潜在 蒸 散发 (BT 从 3 月 中 下 旬 开 


SPEI= 


#2 SPEI 指 数 干 湿 等 级 划分 
Tab. 2 Classification of dry and wet grades of SPEI index 


SPEL <-20 -1.99~-150 -15~-100 -10~-05 -05~05 05~100 10~149 15~200  >2.00 
等 级 。 极端 二 早 WETE PEFR EFEO ”正常 。” 轻 度 湿润 ”中 度 湿润 。 严重 湿润 。 极端 湿润 
(a) 
12 AkKS po 
= OLL 
_ 10} — BST 2 
7 — JLH $ 
2 g] — DHG 5 
a °)— oyza E 
= | 82 | os 
& °)—8 咱 i | 3 
R Ay a R 
= 4 iW | li | E 
* i Wa 机 
Ea is CT a 
= | ERA W 3 
; DIN i y \ NA y 


0 50 100 150 200 250 300 350 
日 时 间 序 列 /d 


0 50 100 150 200 250 300 350 
日 时 间 序 列 /d 


图 2 潜在 蒸 散发 和 实际 蒸 散 发 的 日 变化 


Fig. 2 Daily variation of potential evapotranspiration and actual evapotranspiration 
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台 迅 速 增长 ,6 一 7 月 达到 最 大 ,日 最 大 潜在 蒸 散发 
可 达 11 mm'.d ,随后 逐渐 减 小 ,10 月 以 后 减速 变 
缓 ,总 体 来 看 ,以 7 月 为 中 心 ,对 称 分 布 。 不 同 观测 
点 潜在 蒸 散发 的 日 变化 幅度 有 较 大 差异 ,流域 下 游 
海拔 较 低 的 DHG JLH .CILL 等 地 的 潜在 蒸 散发 变 幅 
很 大 ,海拔 较 高 的 QYZG `QYZ3 、QYZ4 日 变化 幅度 
较 小 。 除 DHG 站 点 外 ,其 余 站 点 的 实际 蒸 散发 
(ET) 的 日 变化 幅度 总 体 上 小 于 洪 在 蒸 散 发 的 变化 
幅度 。 实 际 芋 散 发 中 DHG 观测 点 的 变化 幅度 最 大 ， 
主要 原因 是 该 观测 点 位 于 西 土 沟 流域 下 游 农田 , 受 
灌溉 和 植被 覆盖 度 的 影响 较 大 。 

3.2 流域 年 潜在 蒸 散 发 和 实际 蒸 散 发 量变 化 

由 图 3 可 知 , 西 土 沟 流 域 潜在 蒸 散 发 量 最 高 
达 1478 mm, 最 小 潜在 蒸 散发 量 只 有 426 mm 为 海拔 
4100m 处 的 QYZ4 观 测 点 。 总 体 来 看 , ERE TEAS 
散发 量 小 ,下 游 洁 在 蒸 散发 量 大 。 实 际 蒸 散发 显著 
小 于 潜在 蒸 散发 ,实际 蒸 散发 最 大 的 观测 点 是 下 游 
的 BST 和 机 了 观测 点 ,年 实际 燕 散发 接近 400 mm, 
最 小 的 实际 蒸 散发 量 为 QYZ4 观测 点 处 的 110 
mm。 由 8 个 观测 点 的 潜在 和 实际 蒸 散 发 数据 做 算 
RFE ,得 出 该 流域 年 平均 潜在 蒸 散发 量 为 978 
mm ,年 平均 实际 蒸 散发 量 为 238 mm, 
3.3 不 同 海 拔 梯度 的 潜在 蒸 散 发 和 实际 蒸 散 发 
变化 

8 个 观测 点 涉及 流域 上 、 中 下游, 其 中 流域 下 
游 海 拔 在 1236~1274 m ,平均 海拔 1258 m ,海拔 差异 
小 ,因此 ,将 下 游 的 4 个 站 点 划分 为 同一 海拔 梯度 。 


1750 5 


1500 + e ET, 


1250 4 
1 Vy 


图 3 不 同 观测 站 点 的 年 潜在 蒸 散发 和 实际 蒸 散发 
Fig. 3 Annual potential and actual evapotranspiration at 


different observation sites 
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中 游 阿 克 塞 AKS 海 拔 1686 m, 与 下 游 的 平均 海拔 相 
22428 m, 与 上 游 青岩 子 沟 QYZG 的 海拔 相差 1414 
m, 将 其 与 其 他 站 点 划分 为 同一 海拔 梯度 没有 依据 ， 
因此 ,将 AKS 站 点 单独 划分 为 1 个 海拔 梯度 。 上 游 
的 3 个 站 点 海拔 分 别 是 3100 m、3370 m、4100 m, 将 
这 3 个 站 点 分 别 划分 为 1 个 海拔 梯度 ,因为 3100m 
处 的 QYZG 处 在 山药 位 置 ,3370 海 拔 处 的 QYZ3 站 
点 处 在 较 高 位 置 的 山顶 ,而 4100m 海 拔 处 的 QYZ4 
站 点 位 于 山顶 雪线 以 上 ,此 外 ,这 3 个 站 点 的 下 热 面 
植被 具有 显著 差异 ,因此 将 他 们 单独 划分 进行 分 
析 , 能 够 充分 体现 海拔 对 藻 散 发 的 影响 。 总 结 如 
下 :DHG、BST、CLL、JLH 为 第 1 海拔 梯度 ,平均 海拔 
为 1258 m, AKS 为 第 2 海拔 梯度 ,平均 海拔 为 1686 
m, QYZG 为 第 3 海拔 梯度 ,平均 海拔 为 3100 m, 
QYZ3 为 第 4 海拔 梯度 ,海拔 为 3370 m, QYZ4 为 第 5 
海拔 梯度 ,海拔 为 4100 m。 由 图 4 所 示 TATE ZR 
量 从 低 海 拔 到 高 海拔 逐渐 减 小 ,从 下 游 海 拔 1258 m 
处 到 中 游 海拔 1686 mm 潜在 蒸 散发 下 降 迅 速 ,海拔 
1686 m 到 3370 m 潜 在 蒸 散发 昌 呈 逐渐 下 降 趋 势 ,但 
下 降幅 度 非常 小 ,海拔 3370 m 到 4100m 潜 在 蒸 散 发 
下 降幅 度 较 大 。 实 际 蒸 散 发 量 随 着 海拔 的 升 高 逐 
渐 降 低 ,在 流域 下 游 海拔 为 1258 mm 处 实际 蒸 散 发 最 
高 ,为 356 mm, ,在 流域 上 游 4100 mm 海拔 处 实际 蒸 散 
发 最 低 , 只 有 110 mm。 
3.4 月 尺度 干旱 指数 变化 

计算 西 土 沟 流域 上 、 中 、 下 游 设 置 的 8 个 气象 站 
点 的 潜在 蒸 散发 量 , 根 据 标准 化 降水 蒸 散 指 数 计算 
方法 ,计算 各 个 站 点 的 干旱 指数 SPEI, 分 别 计算 1 个 
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图 4 不 同 海拔 梯度 的 潜在 蒸 散 发 和 实际 蒸 散发 变化 


Fig. 4 Potential and actual evapotranspiration change of 


different elevations 
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月 尺度 和 3 个 月 尺度 的 SPEI 值 ,3 个 月 尺度 的 SPEI 
值 作为 该 年 内 每 个 季度 的 SPEI 指 数 。 由 图 5 可 知 ， 
西 土 沟 流 域 下 游 DHG BST, CLL, JLH 站 点 的 SPEI 
值 从 3 一 9 月 均 为 负 值 ,4 一 9 月 均 出 现 不 同 程度 的 干 
旱 , 其 中 DHG 6 月 和 7 月 出 现 中 度 干 旱 ,BST4 月 和 
9 月 出 现 中 度 干 旱 ,5$ 一 7 月 出 现 极端 干 旱 ,8 月 出 现 
重度 干旱 ;CLL 4 月 出 现 轻 度 干旱 ,5 月 和 8 月 出 现 
中 度 干旱 ,6 月 和 7 月 出 现 重 度 干 旱 ;JLH 3 一 4 月 出 
现 轻 度 干旱 ,5 月 和 9 月 出 现 中 度 干 旱 ,6 一 7 月 出 现 
极端 干旱,8 月 出 现 重度 干旱。 流域 中 游 AKS 站 点 
4 一 9 月 SPEI 指 数 为 负 ,6 月 出 现 中 度 干 旱 ,7 月 为 轻 
度 干旱 。 流 域 上 游 QYZG 和 QYZ3 站 点 分 别 在 5 一 6 
和 8 月 出 现 干旱 ,QYZ4 站 点 全 年 SPEI 为 正 , 较 为 湿 
润 。 图 Si 为 各 个 站 点 的 SPEI 年 均值 ,SPEI 指数 
BST<JLH<CLL<DHG<0<AKS<QYZ3<QYZG<QYZ4。 
3.5 季节 尺度 SPEI 指 数 变 化 

由 图 6 可 知 ,春季 各 个 站 点 的 干旱 指数 均 为 正 
值 且 大 于 0.78, 夏 季 仪 流域 下 游 的 BST、CLL A JLH 
干旱 指数 为 负 , 秋 季 除 QYZ4 站 点 外 ,其 他 站 点 均 为 
负 值 , 且 小 于 -0.3, 其 中 BST 和 JLH 站 点 达到 极端 干 
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旱 , 而 CLL 达 到 重度 干旱 ,DHG 和 AKS 达到 中 度 干 
旱 。 冬 季 BST 和 JLH 站 点 达到 中 度 干旱 ,CLL 站 点 
达到 轻 度 干旱 ,其 他 站 点 属 正 常 或 者 略 湿润 。 综 上 
所 述 ,秋季 是 干旱 发 生 的 主要 季节 ,其 次 是 夏季 。 


4 讨论 


4.1 流域 实际 蒸 散 发 的 估算 及 流域 蒸 散 耗 水 

由 图 7 可 知 ,整体 上 ,上 游 到 下 游 降水 量 逐 渐 减 
小 ,上 游 QYZ4、QYZ3、QYZG 观测 站 点 的 降水 量 平 
均 为 203 mm, 中 游 AKS 站 点 的 降水 为 58.6 mm, F 
游 CLL、BST、JLH、DHG 站 点 的 降水 仅 有 14.6 mmo 
实际 蒸 散发 量 上 游 平 均 为 123.6 mm, 中 游 为 270 
mm, 下游 达 383.7 mm, 上游 到 下 游 实 际 蒸 散 发 量 显 
著 增 大 。 全 流域 具有 上 游 地 区 降水 量 大 于 实际 薰 
散发 ,为 该 流域 的 产 流 区 ,年 降水 量 比 实际 蒸 散发 
量 多 239.8 mm, 上 游 降 水 以 强 降水 为 主 ,降水 强度 
大 ,频次 少 , 因 此 ,降水 容易 形成 洪水 而 流向 中 下 
游 。 中 游 AKS 站 点 蒸 散发 量 比 降 水 量 多 211 mm ,该 
区 下 热 面 为 山 前 洪 积 扇 ,土壤 入 渗 率 大 ,降水 以 及 
上 游 形成 的 洪水 在 该 流域 段 以 及 中 游 末端 的 库 木 
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SPEI 
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图 5 月 尺度 SPEI 指 数 及 各 站 点 年 均 SPEI 值 
Fig.5 Monthly SPEI index and SPEI value of each station 
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DHG BST CLL JLH AKS QYZG QYZ3 QYZ4 
站 点 


图 6 季节 尺度 SPEI 指 数 


Fig.6 SPEIindex change at seasonal scale 


5 降水 
500 ET, 
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观测 点 位 置 


图 7 流域 降水 与 实际 蒸 散 发 的 变化 
Fig.7 Changes of precipitation and actual evapotranspiration 


in river basin 


塔 格 沙漠 地 区 下 渗 ,形成 地 下 径流 。 以 153.05 mm 
为 节点 ,在 流域 上 游 , 当 降 水 量 大 于 153.05 mm 时 , 
实际 蒸 散发 量 低 于 降水 量 ; 在 流域 中 下 游 , 当 降 水 
量 低 于 153.05 mm 时 ,实际 蒸 散发 量 大 于 降水 量 。 
下 游 蒸 散发 量 比 降水 量 多 1476.6 mm , 西 土 沟 流 域 
下 游 位 于 阳 关 镇 ,该 区 域 下 垫 面 复杂 ,有 沙漠 .荒漠 
草地 .荒漠 灌木 .农田 等 ,该 区 的 土壤 蒸发 主要 发 生 
在 沙漠 边缘 .荒漠 草地 等 处 ,植被 蒸腾 主要 发 生 在 
农田 ,灌木 林 区 ,中 游人 渗 的 洪水 在 下 游 流出 地 表 ， 
形成 地 表 径 流 , 供 下 游 农田 灌溉 及 生产 生活 ,此 外 ， 
该 区 的 降水 、 地 下 水 等 均 是 下 游 蔡 散发 的 水 量 来 
源 , 总 体 来 说 ,下 游 是 整个 西 土 沟 流 域 主要 的 水 量 
FERK- 

实际 蒸 散 发 的 观测 比较 困难 ,所 以 ,目前 实际 
蒸 散 发 均 为 估算 值 , 常 见 的 估算 实际 蒸 散 发 的 方法 
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包括 流域 水 量 平衡 法 .Penman-Monteith 公式 - 双 作 
物 和 单 作物 系数 法 .树干 液 流 法 Ts EEA ETE ZS 
渗 仪 法 等 ,由 于 缺乏 实测 资料 ,实际 蒸 散 发 估算 很 
少 涉 及 到 西北 干旱 内 陆 河 流域 ,特别 是 对 荒漠 区 的 
估算 。 本 研究 使 用 在 该 流域 上 中、 下 游 架 设 的 小 
型 自动 气象 站 和 土壤 水 分 观测 设备 实测 的 2018 年 5 
月 1 日 至 2019 年 4 月 30 日 的 气象 和 土壤 水 分 资料 ， 
希望 能 通过 收集 到 的 数据 对 该 区 域 进 行 水 循环 过 
程 中 降水 、 蒸 散发 等 环节 的 定量 研究 ,为 整个 西北 
内 陆 河 浅 山 区 小 流域 的 蒸 散 发 估算 和 和 干旱 特征 研 
究 提供 借鉴 。 本 文 主要 针对 不 同 下 垫 面 .不同 海拔 
开展 流域 内 部 蒸 散发 的 区 域 变化 规律 研究 和 干旱 
特征 的 年 内 不 同时 间 尺 度 的 研究 。 本 研究 使 用 前 
人 关于 作物 系数 K. 的 估算 经 验 公式 “”, 考 虑 到 在 
荡 漠 区 植被 燕 腾 量 很 小 ,主要 是 土壤 蒸发 , 故 在 充 
漠 区 实际 蔡 散 发 佑 算 中 采用 双 作 物 系数 法 ,因为 双 
作物 系数 法 可 以 充分 考虑 到 土壤 蒸发 对 整个 蒸 散 
发 的 影响 ”。 本 人 研究 计算 得 到 西 土 沟 流 域 年 实际 蒸 
散发 为 258 mm, 对 比 已 有 的 关于 干旱 区 内 陆 河 的 实 
际 蒸 散发 研究 , 田 静 等 的 研究 表明 , 近 20 a 来 西北 
地 区 的 年 均 实际 蒸 散 发 为 248 mm, ERX TREE 
玛 纳 斯 河流 域 蒸 散 发 进行 了 研究 ,表明 该 流域 的 年 
实际 蒸 散 发 在 200~250 mm。 李 修 仓 3 的 研究 中 认 
为 ,塔里木 河流 域 的 年 实际 蒸 散发 量 为 278 mm, AS 
文 的 研究 结果 与 他 们 的 研究 结论 较 一 致 ,说 明 Abid 
等 六 和 Rocha 等 下 用 NDVI 来 估算 作物 系数 的 方法 
是 可 以 在 干旱 内 陆 河 流域 应 用 的 。 此 外 , 西 土 沟 流 
域 中 下 游 年 降水 量 不 足 50 mm, 比较 稳定 ,流域 植被 
指数 NDVI 非 常 低 , 中 下 游 的 NDVI 在 0.01~0.47, 平 
均 NDVI 为 0.11, 所 以 ,植被 燕 腾 在 该 流域 所 占 比 重 
显著 小 于 土壤 莱 发 ,土壤 蒸发 的 主要 影响 因素 是 降 
水 和 灌溉 ,对 于 西 土 沟 流 域 这 种 少 降水 、 低 覆 被 的 
地 区 ,土壤 蒸发 受到 限制 ,导致 作物 系数 到 很 小 ,在 
流域 整体 气象 要 素 变 化 较 小 的 情况 下 ,潜在 蒸 散发 
和 实际 蒸 散发 的 年 际 变 幅 较 小 。 结 合 已 有 的 干旱 
内 陆 河 流域 蒸 散发 量 研究 ,本文 计 算 的 实际 蒸 散 发 
量 与 之 较为 接近 ,因此 ,虽然 只 有 一 年 的 气象 和 十 
壤 水 分 数据 ,但 仍然 能 在 很 大 程度 上 代表 流域 乃至 
河西 内 陆 河 浅 山区 小 流域 的 蒸 散发 和 干旱 特征 。 
4.2 干旱 内 陆 河 流域 干旱 指数 SPEI 的 变化 机 制 及 
其 适用 性 评价 

干旱 指数 由 降水 和 潜在 菩 散 发 控制 ,潜在 蒸 散 
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发 主要 由 温度 ,湿度 ,风速 等 要 素 控制 ”六 。 从 上 游 
山区 到 下 游 荒 漠 区 ,降水 温度 .风速 .气压 等 要 素 
发 生 明显 变 化 ,由 此 造成 潜在 蒸 散发 的 区 域 差 异 ， 
从 而 得 出 干旱 指数 的 区 域 差异 。 在 西 土 沟 流 域 ,上 
游 位 于 阿尔 金山 东 缘 ,海拔 较 高 ,温度 气压、 降水 
是 其 与 中 下 游 的 主要 差异 ,低温 和 低压 是 造成 上 游 
山区 潜在 蒸 散 发 量 相 对 较 小 的 主要 原因 。 该 流域 
发 源 于 阿尔 金山 北 项 ,降水 量 大 于 中 下 游 。 而 下 游 
由 于 降水 量 在 50 mm 以 下 ,相对 于 上 游 而 言 , 高 温和 
高 风速 造成 潜在 蒸 散 发 较 大 ,由 此 产生 的 上 游 和 中 
下 游 降 水 与 潜在 蒸 散发 的 差 值 存在 差异 , 上游 差 值 
较 小 ,中 下 游 差 值 较 大 。 根 据 标准 化 降水 蒸 散 指数 
方法 计算 得 出 的 上 中 下 游 的 SPEI 干 旱 指数 存在 差 
异 ( 表 3), 上 游 和 干旱 指数 平均 为 0.38, 中 游 为 0.05, 下 
游 为 -0.53, 根 据 表 2 的 干旱 等 级 划分 ,上 游 和 中 游 
属于 正常 年 份 , 上 游 比 中 游 略 湿润 ,下 游 属 于 轻 度 
于 旱 。 从 海拔 及 下 垫 面 情况 分 析 王 旱 指数 的 变化 
规律 ,海拔 越 高 ,潜在 蒸 散 发 和 实际 蒸 散发 越 低 ”， 
从 上 游 到 下 游 干旱 指数 随 海拔 的 降低 逐渐 减 小 ,上 
中 游 干 旱 指 数 为 正 , 下 游 为 负 。 下 游 SPEI 指 数 降低 
主要 由 湾 在 蒸 散 发 量 控制 。SPEI 指 数 越 低 即 指示 
流域 内 水 量 支出 越 多 ,水 量 收 入 越 少 。 


R3 西 土 沟 流域 上 .中 .下游 年 均 干 旱 指数 
Tab.3 Annual average SPEI index of upper, middle, 
lower reaches of Xitugou watershed 


站 点 位 置 站 点 编号 SPEI 指 数 均值 
上 游 QYZ4 0.86 0.38 
QYZ3 0.07 
QYZG 0.21 
中 游 AKS 0.05 0.05 
下 游 BST -0.84 -0.53 
JLH -0.69 
CLL -0.39 
DHG -0.18 


内 外 采用 的 主流 的 干旱 指数 主要 有 标准 化 降 
水 指数 SPI 帕 尔 默 干旱 指数 PDSI 标准 化 降水 蒸 散 
指数 SPEI, 其 中 PDSI 指 数 在 计算 时 需要 前 期 降水 
量 \ 土 壤 售 水 量 、 径 流 和 潜在 获 散 发 量 等 参数 ,这 些 
参数 有 些 不 易 获 取 , 比 较 依 赖 于 研究 的 区 域 本 身 ”， 
而 SPI 指 数 只 需要 降水 数据 ,没有 考虑 温度 对 干旱 
程度 的 影响 ,SPEI 指 数 是 在 标准 化 降水 指数 SPI 
的 基础 上 加 入 了 温度 对 干旱 的 影响 ,我 国 西 北 干旱 


ChinaXiv 合 作 期 刊 


内 陆 河 流域 ,降水 量 偏 少 ,多 在 50 mm 以 下 ,如果 仅 
考虑 降水 因素 ,会 使 结果 的 准确 度 偏 低 ,因此 ,使 用 
SPEI 指数 能 更 加 合理 的 评估 该 流域 的 干旱 程度 。 
相关 研究 % 也 证 明 ,SPEI 具 有 优 于 PDSI 和 SPI 的 综 
合 优 热 ,而 使 用 Penman-Monteith 公式 计算 参考 作物 
燕 散 发 ,进而 计算 出 的 SPEI 值 (SPEI-PM) 要 比 Harg- 
reaves 和 Thornthwaite 方法 计算 的 潜在 蔡 散 发 计算 
出 的 SPEI 指 数 (SPEI-HA SPEI-TH) 更 能 准确 客观 
AA BCBS NA 。 


5 结论 
(1) 在 年 内 日 时 间 序 列 上 , 西 土 沟 流 域 潜 在 蒸 
散发 和 实际 蒸 散发 以 7 月 为 中 心 , 呈 对 称 变化 。3 
10 月 是 该 流域 蒸 散 耗 水 的 主要 时 期 。 流 域 年 潜在 
蒸 散发 量 为 978 mm, 实 际 蒸 散 发 量 为 238 mms M 
上 游 到 下 游 ,潜在 蒸 散发 和 实际 蒸 散 发 呈现 逐渐 增 
加 的 区 域 变化 规律 。 
(2) 8 个 观测 点 的 年 均 SPEI 指 数 排序 : 拜 水 台 < 
九 莲 湖 < 树 柳 林 < 敦 煌 宫 < 阿 克 塞 < 青 崖 子 3< 青 崖 子 
沟 < 青 岸 子 4, 干 旱 发 生 区 域 : 西 土 沟 流 域 上 .中 、 下 
游 均 会 出 现 不 同 程度 的 干旱 ,从 上 游 到 下 游 , 干 旱 
指数 随 海拔 的 降低 逐渐 减 小 ,上 游 略 湿润 ,到 下 游 
ZWT. PERENE :秋季 是 该 流域 干旱 发 生 
的 主要 季节 ,其 次 是 夏季 ,下 游 的 干旱 程度 和 十 旱 
持续 时 间 要 明显 大 于 中 上 游 。 
(3) 西 土 沟 流 域 上 游 降 水 量 大 于 实际 蒸 散发 
量 ,是 该 流域 主要 的 产 流 区 ,中 游 由 于 下 垫 面 粗糙 ， 
入 渗 率 大 ,地 表 径 流转 化 为 地 下 径流 ,下 游 降水 量 
极 少 ,地 下 径流 出 露地 表 , 供 人 类 生产 生活 以 及 蔡 
散 耗 水 。 流 域 整体 上 莹 散发 量 要 显著 大 于 降水 量 。 
(4) 评 佑 了 利用 NDVI 指 数 计算 的 作物 系数 KK. 
在 河西 内 陆 河流 域 实 际 蒸 散发 估算 中 的 适用 性 , 认 
为 其 计算 出 的 实际 蒸 散发 有 较 好 的 可 信和 度 。 在 SPFI 
指数 计算 中 ,利用 Penman-Monteith 公 式 计算 的 干旱 
此 数 SPEI-PM 更 加 准确 有 效 。 西 土 沟 流域 为 河西 内 
陆 河 浅 山区 的 代表 性 小 流域 ,其 蔡 散 发 规律 及 干旱 
特征 适用 于 整个 无 资料 的 内 陆 河 浅 山 区 小 流域 。 
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Abstract: In order to reveal the evapotranspiration and drought characteristics of typical small watershed in low 
coteau area of Hexi inland river in the arid regions of northwest China, the Xitugou watershed in Dunhuang is used 
as an example to calculate the potential evapotranspiration of eight observation sites by Penman-Monteith formula 
based on the conventional meteorological data. On this basis, the actual evapotranspiration of each observation 
station was calculated by using the method of dual crop coefficient, which was based on NDVI, and then the 
standardized precipitation evapotranspiration index SPEI of the watershed was calculated. The results show that the 
annual potential evapotranspiration and actual evapotranspiration are 978 mm and 258 mm, respectively. From the 
upper reaches to the lower reaches, with the decrease of altitude, the evapotranspiration shows a different degree of 
increase; the duration of drought in the lower reaches of the basin is long, and the whole growth season will have 
different degrees of drought, autumn is the main season of drought; the precipitation in the upper reaches of the 
watershed is greater than the actual evapotranspiration, the actual evapotranspiration of middle and lower reaches is 
greater than precipitation; the K. value of the crop coefficient calculated by the NDVI has good applicability in the 
estimation of actual evapotranspiration, and the drought index SPEI is more advantageous than SPI and PDSI in the 
drought assessment in the arid inland river basin of northwest in China. This study calculated the potential and 
actual evapotranspiration of the Xitugou watershed, evaluated the monthly and seasonal drought characteristics of 
the watershed, and provided guidance for the production and living water consumption of the watershed, especially 
in the middle and lower reaches. This research obtained the general disciplines to no data river in the low coteau 
area of the whole Hexi inland river. 
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